From adiabatic piston to non-equilibrium hydrodynamics 
Ken Sekimoto 1,2 , Antoine Fruleux 2 ' 1 , Ryoichi Kawai 3 and Nathan Ridling 3 

'Matieres et Systemes Complexes, CNRS-UMR7057, Universite Paris-Diderot, 75205 Paris, France, 
2 Gulliver, CNRS-UMR7083, ESPCI, 75231 Paris, France, 
3 Department of Physics, University of Alabama at Birmingham, 



Based on the new concept of the momentum transfer deficiency due to dissipation (MDD), the 
physical basis of the mechanism of "adiabatic piston" is explained. The implication of MDD in 
terms of hydrodynamics under non-equilibrium steady state also discussed. 
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1. Introduction 

Since the 1920's a simple question associated to non-equilibrium statistical physics has been 
addressed: If we put a Brownian piston of mass M between the two semi-infinite cylinders each be- 
ing filled with an ideal gas consisting of particles with mass m(<sc M), what is the non-equilibrium 
steady state ? Here the temperature and pressure of the gas in the left cylinder are prepared at (T, p) 
while those in the right cylinder are at (T' ,p). The surface area of the piston is the same on both 
sides. It is clear that, if the piston were firmly fixed and if the piston is "adiabatic", then there would 
be no net force on the piston because the gas in each cylinder remains in equilibrium and presses 
the piston by the same pressure but in the opposing directions. When the Brownian motion of this 
adiabatic piston is allowed, however, this motion will allow the energy transfer from the hotter gas 



(e.g. the side of T if T > T') to the colder gas (ibid. T') [Feynman et al.(1963)(§46.1-§46.9) |. 



The question is if this non-equilibrium precess leads to a non-vanishing net force on the piston. 



The macroscopic thermodynamics cannot answer this type of question [Callen(1965)|, neither the 



Langevin description can give an answer to this type of setup [Van den Broeck et al.(2004)| . While 
the stochastic energetics | |Sekimoto(1997)| |Sekimoto(2010)] can describe correctly the heat flow, 
the non-equilibrium force is beyond the resolution of this level of description. 

Many calculative studies have been reported in the past both on this problem and also on a 
class of Brownian ratchet models, which turned out to be essentially the same problem as adi- 
abatic piston [ Fruleux et al.(2012)] . All these studies have been done using either by ad hoc 



treatment of Master-Boltzmann equations with truncated moment hierarchical expansion with e 
■^m/M as small parameter, or by molecular dynamics (MD) simulations, see the references cited 
in [Fruleux et al.(2012)j. Both the perturbative approach and the MD simulation consistently con- 



cluded that the Brownian piston will move steadily towards the hotter gas. Nevertheless a clear 
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physical understanding was missing. A frequently given hand-waving argument was that the hotter 
side losing the heat has locally lower pressure. But it is not a valid argument. If the ideal gas 
is used, the cooled particle will never hit again the Brownian piston while the freshly colliding 
particles are characterized by the equilibrium parameters, (T, p) or (T',p). The adiabatic piston 
has, therefore, remained among "Some problems in statistical mechanics that I would like to see 
solved" |Lieb(1999)[ . 

It is only very recently | |Fruleux et al.(2012)jj that the physical explanation to the above ques- 
tion was definitely given. The key is to take into account the inteiplay between the energy transfer 
and the momentum transfer at the gas-piston interfaces. Once this point is understood, the results 
of elaborated perturbative calculations could be perfectly reproduced just by a few lines' calcula- 
tions, except for an overall numerical factor. The purpose of the present paper is to summary the 
basic idea and discuss its generalization. The organization of the paper is the following: In the next 
section (§ [2]) we recapitulate the main line of this mechanism. Especially the key concept of the 
momentum transfer deficit due to dissipation (MDD) is explained using a simple argument. The 
relation to the traditional calculative approach is also mentioned. As a prologue to the extension 
to the dense gas case, we introduce in § 13.21 a toy model that shows how the energy flow and mo- 
mentum flow having different symmetry in space and time can make the uniform pressure and the 
heat conduction compatible. In §[33] we show the implication of the MDD to the non-equilibrium 
hydrodynamics of dense hard-core gas. 



2. Review of the physics of adiabatic piston 

We outline the concept of MDD as the underlining mechanism of the adiabatic piston. The 
readers might refer to [ Fruleux et al.(2012)] [ Kawai et al.(2012) | for the technical details and its 



generalization to inelastic case. 

The essential point of the adiabatic piston is more clearly grasped when the Brownian piston 
separating the ideal gases is trapped by a potential force such as an elastic spring (Fig. \Q so that the 
mean velocity of the piston vanishes ({V) = 0). If there appears a non-equilibrium force Fness on 
this trapped Brownian piston, the steady state velocity (V) of the piston in the absence of trapping 
is given by the balance with the passive frictional force, Fness _ (y + y')(V) = 0, where y and y' 
are the friction coefficient of the Brownian piston against the respective gas. 




Fig. 1 . Trapped adiabatic piston. 



The first step is to realize that the Brownian motion of the piston serves merely as the medi- 
ator of the energy transfer, or heat, from hotter gas to the cooler gas. While correlation between 
collisions with the piston by the hot gas particles and by the cold gas particles are essential for the 
irreversibility of this purely mechanical problem, we can bypass all the details for the purpose of 
understanding the non-equilibrium force Fness- The rate of the energy transfer per unit surface of 
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the piston, jr> can be found using the stochastic energetics ||Sekimoto(2010)l or even by a heuristic 



argument | Parrondo and Espanol(1996)] . The result reads, / e) = (k#T _ k^T')/[M(y 1 + y' *)], 



where y = y/A and y' = y' /A with A being the area of each piston surface. 

Once we know the energy transfer rate across the gas-piston interface, we can concentrate on 
the following problem: When an ideal gas prepared in the equilibrium characterized by (T, p) is 
brought into contact with the wall that absorbs (or injects) energy at the rate / e) per unit surface 
(Fig. |2]), how the pressure on the wall is modified from pi Consider that gas particles with a 



P(0) 



EeqC*), P 



Fig. 2. Gas in contact with an energy-transferring wall at x — 0. The contact density is denoted by p(0). In 
general the effective temperature r new (x) should depend on the position x. 

typical velocity component normal to the wall v± collides the energy-transferring wall. They are 
reflected back with a velocity v^ ut . While the typical incoming velocity should be the thermal 
velocity = v tn = yfkoTjm, except for the numerical factor, the typical outgoing velocity depends 
on the energy transfer rate ft e > and the collision rate v co j per unit area on the energy-transferring 
wall through the energy balance condition: 



Vcol 



_ m 2 

" 2 V in 
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- jv out 



Assuming that the energy transfer is sufficient small, the right hand side (r.h.s.) of the above 
equation is approximated as 



v ± + W 1 - A 

i x i _i_ i \ in 1 out i , x i _l i\ 
(mv- - m|v out |) - (mv- - m|v out |)v th . 



Substituting this result into the above equation, we find 

j(e) 

(mv± - m|v^ ut |)v co i - — 

Vth 



(1) 



The left hand side (l.h.s.) of this relation gives the momentum transfer deficit due to dissipation 
(MDD). In other words, upon the collision, the gas particles kicks the wall less strongly in non- 
equilibrium than equilibrium if a part of their incoming kinetic energy was taken out by the energy- 
transferring wall. In terms of the net momentum transfer rate per unit surface, j ( P ) - LJ - - v co i(mvth + 
l v outD' tne equilibrium value, p = 2rav t hV co i, is corrected by this MDD to give 

/(e) 
Vth 
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In retrospect, the traditional approach through the Master-Boltzmann equation could have given 
the same insight. For the setup of Fig. [2] this equation can be written as follows: 

d t P(X, V, t) = -Vd x P(X, V, t) - [-y fc V - d x U(X)]d v P(X, V, t) 

- f W(V'\V)P(X, V,t)+ f W(V\V')P(X,V' ,t) + ^^-d\P(X,V,t), 

JV Jv 7b 

where P(X, V, t) is the probability density of the position X and velocity V of the wall as a Brownian 
piston, and U(X) represents the trapping potential energy. The heat absorption by the wall is 
modeled by the coupling to a Langevin bath | |Sekimoto(2010)l at the temperature Tb with the 
coupling, i.e. friction, constant ji,. The collision of the gas particles is represented by the velocity 
transition rate, W(V given by 

VYYl m I vyi + M\ 

IkkbT \ 2m l 

where A is the surface area of the wall, v 1 is the normal component of the incoming velocity of gas 
particle, and H(z) is the Heaviside unit step function, v 1 is the function of wall's velocities before 
(V) and after (V) the collision, respectively, through the momentum conservation rule, 

2m , 

V = V + (v x - V). 

m + M 

The truncated equations for the first two moments of V can be derived from the above Master- 
Boltzmann equation, and the results read 

d(V) I i (e) \ 

M-y- = -(U'(X)) - (y + y b ){V) + [p-c J —]A 
dt \ vth / 

M d{v2 = _ {VU > {X)) _ z [<y 2 > _ kT] _ zw } _ kT] + c , (v) 

dt MM 



where c = and the other constant c' is irrelevant as we shall see immediately below. In 

the steady state, not only (V) = d/dt = but also (VU'(X)) vanishes. Then the second moment 
equation tells that the kinetic temperature of the Brownian piston, = M(V 2 ), is given by the 
well known formula of Langevin dynamics, 

; r _ 7 k ^ T + 7b^T b 

KB-* kin - 

7 + 7b 

Moreover, the second and the third terms on the r.h.s. of the second moment equation gives the 
energy transfer to [from] the wall, respectively: 

j (e) = -JjW 2 ) ~ k ^T] = jjl(V 2 ) - k B T b ]. 
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With p e > thus known, the first moment equation in the steady state is nothing but the momentum 
balance condition: _ ^ 

-(U\X)) + \p-c^-jA = 0. 

Our physical reasoning, therefore, reproduces completely the traditional result except for the nu- 
merical factor c. Moreover, our explanation allows to treat the adiabatic piston, Brownian ratchet 
models [Van den Broeck et al.(2004)| , or inelastic piston | |Costantini et al.(2008)| ] on the same foot- 
ing lFruleuxetal.(2012)l. 



3. Momentum transfer of a gas with heat transport 

3.1. Preliminary argument 

The mean free path £ map of an ideal gas is infinite because the particles undergo no collisions. 
Knudsen number Kn = £ ma p/L sys is therefore infinite with any system size, L sys . The macroscopic 
thermo-hydrodynamics | |Landau and Lifshitz(2000) | supposes the opposite limit, Kn «; 1 . When 



we study the thermo-hydrodynamics with energy-transferring boundaries, the physical ideas ob- 
tained in the previous section should, therefore, be applicable only to the vicinities of those walls 
probably with some modifications. The main question is how to reconcile the formula Eq. [T]for the 
ideal gas with the macroscopic description of thermo-hydrodynamics with non-equilibrium bound- 
ary condition. In this paper we limit ourselves to the steady states with vanishing macroscopic 
velocity of the gas. The conservation laws of mass, momentum and energy then impose the con- 
stancy of those fluxes. Below we study first by a purely mechanical toy model that shows the basic 
compatibility between these fluxes and their nature of symmetry in space and in time (§§E2j). Then 
we go onto the dense hard-core gas with Kn «; 1 (§ 



3.2. Toy model 

We begin by a very elementary kinetic model to discuss the interplay of the energy and momen- 
tum transfeiQ. We take up a single gas particle on the jc-axis bounded by the energy-transferring 
walls at x - and at x = L sys , which are macroscopically fixed in space, see Fig. [3] We further 
simplify that the hot wall (on the left) receives the particle of velocity -v_ and returns with the 
velocity v + with < v_ < v + . The cold wall (on the right) does the opposite operation. The micro- 
scopic mechanism underlying these reflections are irrelevant for our argument. (On might imagine 
the two tennis player engaging in a rally.) 

Before counting the momentum and energy flux, / e) and f p \ we impose the vanishing of the 
mass flux, f m) in the steady state: 

J - (p+mv + - p-mv-)x = 0. 
This is satisfied by the densities of rightward and leftward particle, respectively, 

v+ 

(v+ + v„)L sys 

1 This is a simplified version of Knudssen heat transfer, see, for example, jStruchtrup(2005) ], page 25. 
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time 




Fig. 3. Space (x)- time trajectory of a particle between the hot (x = 0) and cold (x = L sys ) walls 



Also the collision frequency on each wall, v co i is found to be 

V+V- 



Vcol 



With this Vcoi, the energy transfer rate / e) is 



(v+ + v_)L 



sys 



m 



v+v_ 



m 

(v + +v_)L sys 2 



v-)- 



^sys 



while the momentum flux reads 



V+V- 



(v+ + v_)L s 



sys 



— JUT. 

-^sys 



We verify that is odd under time or space-inversion, while / p) is even under these operation. 
The symmetry of / p) can also be seen from Fig [3] This model, although simple, shows how the 
directed energy transfer is established without gradient of momentum flux. In other words, the 
pressure on the hot and cold walls are the same. 

To see more in detail the process at the walls, we refer to the contact value theorem, p = 
k^TpiO) [Henderson et al.(1979)[ , which give the equilibrium momentum transfer to a hard wall 
by a hard core gas in terms of the equilibrium temperature T and the gas particle density at the 
closest contact surface of the hard wall. Our interest is the case with energy-transferring walls, see 
Fig. 12 When the walls transfer the energy, the p should indicate the total momentum flux f p \ i.e., 

Aieq = \] {P) \ = m(v+ + V_)v col 

and k^T ne nl2 the kinetic energy per particle, 



^B^neq m p+V + 2 + p_V_ 2 



m 

= — V+V- 



2 2 p + + p_ 2 
Since the total density p(0) on the wall is, by the homogeneity, 

/ 1 1 \ 1 

p(0) = — + — y col = — 
\v+ v-j L sys 
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we arrive at a form of the contact value theorem in non-equilibrium. 

Pneq = k B T Deq p(0). (2) 



In other words, while the symmetry allows the correction to the r.h.s. of the form oc (v + - v_) or 

2 2 

oc js e > , the contact value theorem holds up to the order of 0(j {e) ) if p eq and T eq are appropriately 



chosen (cf.|Komatsu et al.(2008) 1). 



3.3. Non-equilibrium hydrodynamics 

In the non-equilibrium steady state with heat flux of a dense hard-core gas with Kn «; 1, the en- 
ergy flux vector field, f e \ and momentum flux tensor field, j (p) , must satisfy the basic conservation 
laws: 

v • f e) = o, v • f p) = o. 

If the wall is perpendicular to the ^-axis, the system is homogeneous in y and z directions and the 
above conditions are reduces to 

= const. \f x = const. 

If the heat conduction obeys approximately the Fourier's law, = kjVT, with kj being the heat 
conductivity, the temperature gradient keeps constancy of the energy flux. As for the momentum 
flux f p \ the symmetry argument or Curie principle [ Curie(1894)] allows the anisotropy of the 
type = pi + a(xx - il) with a characterizing the deviatoric part of the flux due to heat flux 

||jc. However, seeing that the (9(/ e ^ 2 ) contribution was missing in the above simple model Eq. |2j 
we simply identify = p neq l to be the pressure in the present approximation. If p neq obeys 
approximately the equilibrium equation of state, p = p{p, T), among the pressure p, temperature 
T and the density p, the density p(x) varies in a manner locally compensating the heterogeneity of 
the temperature T(x) so that the p neq remains homogeneous. 

Our concern is how we can relate the momentum flux and the energy flux in the dense hard- 
core gas where p neq reflects already both the incoming and outgoing particles. Below we will 
indicate that the relation like Eq. Q] corresponds to the skewness of the velocity distribution of 
particles (especially) at the energy-transferring wall. To be concrete we imagine the dense hard- 
core gas which is conducting the heat rightwards up to the energy-transferring wall at x = without 
convection (Fig. [2]). We also assume that the wall exchanges only the ^-component of momentum. 
Now we introduce the velocity distribution function f(y x ; x) per unit volume of gas particles. Then 
the particle density p(x) is given by 



p(x) = J f(y x ;x)dv x . 

The conditions on the fluxes of mass, momentum and energy along the x axis are given, respec- 
tively, as 
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Pneq = * ■ \ (p \x) ■ x = J mv 2 J(v x ;x)dv x 
= f e \x)-x = J jv 3 x f(v x ;x)dv x , 



(3) 



where p neq and are independent of the position x. 

Now we focus on the thin slab of the distance «; ( m f p from the energy-transferring wall. In this 
slab we assume that the gas particles undergo practically no collisions except for with the wall. We 
introduce the partial momentum fluxes associated to the incoming particles, j.^ and to the outgoing 
particles, f^ v right before the wall: 



J 



.v=0- 



= f mv 2 x 



l : H(+v x )f(v x ;x)dv x , 



M 

" out 



X=Q- 



= ^ mv 2 x H(-v x 



)f(v x ;x)dv x . 



In equilibrium where f(v x ;x) are symmetric with respect to v x , the both partial fluxes are the same. 
In the presence of the heat flux it is no more the case. While the asymmetric velocity distribution 
for Knudsen heat transfer, i.e. the above toy model, is usually singular and far from Maxwellian, 
the collisions make the velocity distribution looks like skewed Maxwell distribution. We, therefore, 
assume an approximate forrro 



f(v x ;x)\ x ._ 



+ C 3 



CT 



For a week energy flux, the terms containing c\, c% and C3 are regarded to be small perturbations 
with respect to the main term cq. The expression of the density and the flux conditions mentioned 
above impose 



Pi 



r=0- 



V2^ 



Pneq 

V27TW? 

2 fe) 



(co + C 2 )cr 
C\ + 3C3 

(c + 3c 2 )o- 3 



(ci + 5c 3 )cr 5 = 2c 3 cr 5 , 



(4) 



3 V27T»J 

where in the last line we used the vanishing mass flux condition in the second line; c\ - -3c 3 
Finally the difference between the partial momentum fluxes j. and read 



2m 



= (ci + 4c 3 )cr 4 = c 3 cr 4 , 



where again we used the vanishing mass flux condition. Then if we introduce the squared average 
of the particle velocity (noting \c2\lcQ «: 1 for weak non-equilibrium), 



vilx=o- 



T^neq _ 1 + 3c2/cp 
PUO- 1 + C 2/C0 



-cr 



cr 



2 On the wall, x = 0, the very MDD implies the discontinuity in f(v x ;x) at v r = 0. Here, however, we shall use a 
smoothed form as qualitative model. See also, for example, [Struchtrup(2005)J, page 202. 
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we obtain the relation of MDD and the energy flux, reminiscent of Eq. [T] 



(&>-&)\ = C J - 



V x \x=0- 



1/2 



with a numerical factor C = V2/(3-0r)> which is subject to our approximations. 

In onclusion, the dense hard-core gas conducting the heat carries also momentum through the 
asymmetric distribution of particle's velocity, and the concept of MDD is a neat way to explain the 
relation between the energy flux and the partial momentum fluxes in the non-equilibrium steady 
state. 



We would like to acknowledge the organizers of the 25th Smoluchowski Symposium. 
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